arXiv:1503.05965v2 [math.CA] 29Jul2015 


Calculus in the ring of Fermat reals 

Part I: Integral calculus 

Paolo Giordano^’Enxin 

University of Vienna, Austria 
^ University of Vienna, Austria 


Abstract 

We develop the integral calculus for quasi-standard smooth functions defined 
on the ring of Fermat reals. The approach is by proving the existence and 
uniqueness of primitives. Besides the classical integral formulas, we show the 
flexibility of the Cartesian closed framework of Fermat spaces to deal with 
infinite dimensional integral operators. The total order relation between 
scalars permits to prove several classical order properties of these integrals 
and to study multiple integrals on Peano-Jordan-like integration domains. 
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1. Introduction 

In spite of old controversies concerning the existence and rigor of in¬ 
finitesimals, it is remarkable that using only elementary calculus, one can 
define and deeply study non-trivial rings containing the real field and in¬ 
finitesimals. This goal has been achieved e.g. with the Levi-Civita field 
([21, 22, 2, 3, 25, 26]), the Colombeau ring of generalized numbers ([4, 5, 24]) 
or the ring of Fermat reals ([9, 10, 12, 14, 8]). 

A common erroneous opinion is that in non-Archimedean theories, we 
pay the price of needing a non-trivial knowledge of mathematical logic ([16]) 
or switching to intuitionistic logic ([18]) or using some non-trivial set theo¬ 
retical methods ( [6] ). The non-Archimedean theories mentioned above need 
none of these. Another common feeling is that, to deal with infinitesimals 
one needs to have a great formal attitude in doing mathematics. Of course, 
there are half-way approaches trying to require a less formal path: e.g. at 
least up to a certain level, also nonstandard analysis, surely the most power¬ 
ful theory of infinitesimals and infinities, can be presented without the need 
to face formal logic ([23]). 

Aside from theories based on analysis, several algebraic methods are also 
available ([27, 2, 19, 1]). However, in those cases, the interested reader has to 
be able either to proceed formally, or to develop an intuition related to such 
algebraic infinitesimals. For example, in [10] it is argued that an intuitive 
picture of infinitesimal segments of lengths h and k such that = 0 

but h-k 0 is not possible. These numbers are frequently used in Synthetic 
Differential Geometry (SDG, [18]). At the end, it is a matter of taste about 
what approaches are felt as beautiful, manageable and in accordance with 
our philosophical approach to mathematics. 

The ring *M of Fermat reals, can be defined and studied using only el¬ 
ementary calculus; see ([9]). It extends the field M of real numbers, and 
contains nilpotent infinitesimals, i.e., h G *M^o such that = 0 for some 
n G N>i. The methodological thread followed in the development of the 
theory of Fermat reals has been always the necessity to obtain a good dia¬ 
logue between formal properties and their informal interpretations. Indeed, 
to cite some results of this dialogue, we can say that the ring 'M is totally 
ordered and geometrically representable ([10]). 

Related to the ring ’M, we have a good notion of smooth function called 
quasi-standard smooth function. These functions are arrows of the category 
*C°° of Fermat spaces ([12, 8]), which is essentially a generalization of the 
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notion of diffeological space ([17]), but starting from quasi-standard smooth 
functions defined on arbitrary subsets S C instead of smooth functions 

defined on open sets t/ C The category is a quasi-topos [15], and 

hence has several desired categorical properties, such as Cartesian closed¬ 
ness. Furthermore, *C°° contains finite-dimensional smooth manifolds, and 
inhnite dimensional spaces like Man(M, A^), the space of all smooth func¬ 
tions between two smooth manifolds M and N, and integral and differ¬ 
ential operators, etc. One of our goals is to develop differential geometry 
on these spaces using nilpotent infinitesimals. For the first steps, we need 
to develop differential and integral calculus, which is also preliminary to 
other advanced topics, such as calculus of variations, generalized functions, 
stochastic infinitesimals, etc. 

As a start, in the present work, we develop integral calculus for quasi¬ 
standard smooth functions. Later in another paper, we will develop differ¬ 
ential calculus. This exchange is motivated by the need of using Taylor’s 
theorem with integral remainder in differential calculus. 

In our approach, the derivatives are defined using the existence and 
uniqueness of the quasi-standard incremental ratio (see the Fermat-Reyes 
Thm. 2), and the integrals are defined using the existence and uniqueness 
of primitives (see Thm. 12 below). Following this approach, we have an 
extension of the classical results of calculus and also a quick agreement with 
several informal calculations used in physics (see e.g. [10]). 

What are the new results following this approach, compared to classical 
analysis^ or NS A? Besides the possibility to use nilpotent infinitesimals to 
simplify several calculations, we will see that the Cartesian closed framework 
of Fermat spaces permits to easily prove that all the classical infinite 
dimensional integral and differential operators are smooth maps, i.e., they 
are morphisms in *C°°. On the other hand, convenient vector spaces and 
related smooth functions ([7, 19]) can be embedded both in the category Dig 
of diffeological spaces and in the category of Fermat spaces (see [12, 17]). 
Therefore, one of the best calculus in locally convex topological vector spaces 
is still available in *C°°. Finally, the Fermat functor *(—) : Dig —>■ 

([8, 15]) connects these categories and has very good preservation properties. 

What are the new results compared to smooth infinitesimal analysis 
(SIA; see e.g. [20])? In the ring of Fermat reals, we have a total order rela¬ 
tion rather than only a pre-order, and indeed we are able to prove classical 
inequalities and order properties of integrals. Because of the simplicity of 
the whole model (if compared to those of SIA), in our framework we can 
not only formally repeat several proofs of SIA, but also suitably generalize 


®Here by classical analysis, we mean the most developed theory of infinite dimensional 
spaces, i.e., locally convex topological vector space theory. We recall that the category of 
these spaces is not Cartesian closed; see [11]. 
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some classical existence theorems like the inverse function theorem and the 
intermediate value theorem. 

The structure of the paper is as follows. 

In Section 2, we summarize the basics of Fermat reals: unique decompo¬ 
sition of elements, Taylor’s formula with nilpotent increments, Fermat topol¬ 
ogy, extension of ordinary smooth functions, well-ordering, quasi-standard 
smooth functions, etc. 

In Section 3, we review another topology on Fermat reals called the 
cu-topology ([14]), which is only used in this section. Unlike the Fermat 
topology, we show that in general, quasi-standard smooth functions are not 
continuous with respect to this topology. On the other hand, we prove a 
uniqueness theorem (Thm. 5), which says that every quasi-standard smooth 
function is determined by its values on any dense subset of the domain in the 
cu-topology. This uniqueness theorem is a generalization of previous results 
like the cancellation law of non-infinitesimal functions ([12, Lem. 25]), and 
it will be useful for subsequent sections. 

In Section 4, we show that every interval and its interior in "M is con¬ 
nected (Lem. 6), and we use this to prove the constant function theorem 
(Thm. 7), i.e., that only constant quasi-standard smooth functions have zero 
derivatives on a non-infinitesimal interval (see Sec. 2). From now on, we 
require the ambient quasi-standard smooth function to be defined on a non¬ 
infinitesimal interval, due to the uniqueness of the incremental ratio in the 
Fermat-Reyes theorem (Thm. 2). Since the domain of every quasi-standard 
smooth function defined on an interval of *M can always be extended to a 
Fermat open set (although the extension may not be unique; see Lem. 8), 
we generalize the Fermat-Reyes theorem for quasi-standard smooth func¬ 
tions defined on a non-infinitesimal interval instead of a Fermat open set, 
and hence introduce the notion of left and right derivatives (Lem. 10). 

In Section 5, we prove the existence and uniqueness of primitives for 
quasi-standard smooth functions (Thm. 12 and Thm. 13) by the sheaf 
property of such functions and the constant function theorem (Thm. 7). 
From this, we define integrals of quasi-standard smooth functions (Def. 14), 
which is a natural generalization of integrals of ordinary smooth functions 
(Rem. 15). Note that, although we require the ambient quasi-standard 
smooth function to be defined on a non-infinitesimal interval, the integrals 
can be defined on an inhnitesimal interval. From the proof of Thm. 12, we 
get a useful expression of integrals (Lem. 16) as a sum of integrals of the 
local expression of quasi-standard smooth functions as extension of ordinary 
smooth functions with parameters. 

Although most results about integrals of quasi-standard smooth func¬ 
tions can be proved using Lem. 16, the local expression, (see Rem. 18), we 
prefer “global” proofs. In Section 6 and the first part of Section 7, we prove 
the classical formulas of integrals calculus (Thm. 17, Thm. 19, Lem. 23 and 
Cor. 24). The “global” proof of the additivity of integrals with respect to 
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the integration intervals (Cor. 24) relies on the commutativity of differenti¬ 
ation and integration (Lem. 22), and the latter is a consequence of the 
smoothness of some infinite dimensional integral operators (Lem. 20 and 
Cor. 21). Finally, we get a form of mean value theorem called Hadamard’s 
lemma (Cor. 25). 

In the rest of Section 7, we study the standard and infinitesimal parts of 
an integral (Def. 26 and Lem. 27), and we review divergence and curl using 
suitable (nilpotent) infinitesimals in Fermat reals. As a result, any integral 
over an inhnitesimal integration interval is infinitesimal, and any integral of 
an inhnitesimal-valued quasi-standard smooth function is infinitesimal. 

Thanks to the total order on 'M, in Section 8, we get some similar 
inequalities for our integrals as the classical ones, with respect to arbitrary 
integration interval (monotonicity of integral (Lem. 28 and Thm. 29), a 
replacement of the inequality for absolute value of quasi-standard smooth 
functions (Thm. 30), and the Cauchy-Schwarz inequality (Thm. 31)). 

In Section 9, we study multiple integrals of quasi-standard smooth func¬ 
tions, with integration domain a finite family of pairwise disjoint boxes 
(Thm. 34 and Thm. 37), and we prove Fubini’s theorem (Thm. 36). 

2. Background for Fermat reals 

In this section, we summarize the basics of Fermat reals by listing some 
properties that permit to characterize the ring *M (see [13] for a formal 
list of axioms). For a more detailed presentation and for a description of 
the very simple model, see [9, 8]. See also [14] for the study of *M using 
metrics, the characterization of its ideals, roots of nilpotent infinitesimals, 
some applications to fractional derivatives, and a computer implementation. 

At the end of this section, we hll in a gap in the proof of the Fermat- 
Reyes theorem in [12], so that the statement of the theorem (see Thm. 2) is 
slightly different, but the differential calculus developed in [12] still holds. 

Fermat reals can be defined as a quotient ring of the ring of little-oh 
polynomials. Equivalently, every Fermat real x G *M can be written, in a 
unique way, as 


N 



( 1 ) 


2 = 1 


where °x, ai, G M are standard reals, ai > 02 > ■ ■ ■ > otv > 1, Oi 7 ^ 0, 
and dfa verifies the following properties: 


dta ■ dtb = df ab 


( dta)^ = dta Vp G M>i 

V — 

dta = 0 Va G M<i. 
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The expression (1) is called the decomposition of x, and the real number 
°x its standard part. The greatest number ai =: uj{x) is called the order of 
X and represents the greatest infinitesimal appearing in its decomposition. 
When X G M, i.e., x = °x, we set a;(x) := 0. We will also use the notations 
uJi{x) := ai and °Xi := ai for the i-th order and the i-th standard part of x. 
The order a;(—) has the following expected properties: 


uj{x + y) < max{a;(x),a;(?/)} 
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H-^ if °x = °y = 0 and x ■ y ^ 0. 


uj{x-y) uj{x) uj{y) 

Using the decomposition of x, it is not hard to prove that for k G N>i, 
x^ = 0 iff w(x) < k. For k G M>o U {oo}, the ideal 

Dk := {x G *M I °x = 0, a;(x) < A: + 1} 


plays a fundamental role in k-th order Taylor’s formula with nilpotent in¬ 
crements (so that the remainder is zero). Indeed, for k G N>i we have 
that Dk = {x G *M I = 0|. We simply write D for Di. Every ordinary 
smooth function / : A — > M, dehned on an open set A of can be 
extended to a function ' f ■.* A —)■ 'M defined on 

•A := |x G •M‘^|°x G a} , 

preserving old values /(x) G M for x G A. Note that, for A,B open sets 
of Euclidean spaces, we can naturally identify *(A x B) with *A x *B, and 
hence there is no ambiguity to write for (‘M)'’* or see [12, Thm. 

19]. The mentioned Taylor’s formula is therefore 

VfteDf: •/(x + /.)= (2) 

\j\<k 


where x G and Df = x .. . xDj.. 

There is a natural topology on "M'’* consisting of sets of the form * A for 
open sets ^4 C M^. We call this topology the Fermat topology, and open sets 
in the Eermat topology the Fermat open sets. Unless we specify otherwise, 
in the present work we always equip every subset of with the Eermat 
topology when viewed as a topological space. 

It may seem difficult to work in a ring with zero divisors. But the 
following properties permit to deal effectively with products of nilpotent 
inhnitesimals (typically appearing in Taylor’s formula of several variables) 
and with cancellation law: 


n 
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X is invertible °x ^ 0 

(If X • r = X • s in *M, where r, s € M and x 7 ^ 0) r = s. 

The (commutative unital) ring *M is totally ordered, and the order rela¬ 
tion can be effectively decided starting from the decompositions of elements 
as follows. For x € *M, if °x 7 ^ 0, then x > 0 °x > 0; otherwise, if 

°x = 0 , then x > 0 °xi > 0 . 

For example, dt^ — 3 dt > dt > 0, and dt^ < d4 if a < 6 . In the 
following, we will use symbols like [a,b] := {x G *M | a < x < 6 } for intervals 
in *M, whereas intervals in M will be denoted by symbols like [a, ft]® := 
[a, 6]nM. By an interval in ’M we mean a set of the form [a, b], (a, 6), (a, b] or 
[a, 6 ), where a <b. An interval is called non-infinitesimal if °a < °b. Since < 
is a total order on ‘M, it is easy to prove e.g. that a = inf (a, b), b = sup(a, b). 
Therefore, each inteval of *M uniquely determines its endpoints. 

The calculus we begin to develop in the present work concerns quasi¬ 
standard smooth functions, which were introduced in [12]. Recall that a 
function / ; S —^ T with S C and T C •R'', is called quasi-standard 
smooth if locally (in the Fermat topology) it can be written as 

f{x) = *a{p,x) Vx € n 5, (3) 

where a G C°°{P x R, M'-) is an ordinary smooth function defined on an 
open set P x V C RP x R^, and p G *P is some fixed p-dimensional pa¬ 
rameter^. Every quasi-standard smooth function is continuous when both 
domain and codomain are equipped with the Fermat topology (see [12, Thm. 
13]). We write S*R°°(A,R) for the set of all quasi-standard smooth func¬ 
tions from A C ’R^ to R C •R*’. Moreover, when A and B are open sets 
of Euclidean spaces, S*R°°(A, R) = ‘C°°{A,B) = C°°{A,B) (see [12, Thm. 
18(4)]). Therefore, the notions of “ordinary smooth function” (i.e., func¬ 
tions such that all iterated partial derivatives exist and are continuous), of 
“quasi-standard smooth function” and of “arrow of the category (see 

[12, Def. 16]) are consistent. From now on, we call every such function 
a smooth function, and one can tell which kind it is from its domain and 
codomain. 

In [12], the Fermat-Reyes theorem, which is essential for the development 
of the differential calculus on *R, was presented. Unfortunately, there is a 
gap in formula ( 22 ) of [ 12 ], since nothing guarantees that the neighborhood 
*R of h is sufficiently small for that formula to hold. Nevertheless, it is 
clear that we need to guarantee the validity of the Fermat-Reyes formula 
(see ( 6 ) below) only for h sufficiently small. In the following definition, we 
precise the meaning of the sentence “for h sufficiently small, the property 
V{h) holds”. 


^Note the different fonts used e.g. for the point p £ P and the dimension of R'’ 3 P. 
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Definition 1. Let V^h) be a property of /i € 5 C ‘M. We write 

: V{h), 

and we read it as for h & S sufficiently small V{h) holds, if 
3p € ]R>o V/i G 5 : |/i| < p P(/i). 

Note explicitly that this notation also includes the special case 5 C M. Note 
also that, in a formula like Vx V^/i : 'P{x, h), if we say that V{x, h) holds for 
\h\ < p, then p depends on x. 

For example, if [/ C is a Fermat open set and v G *M'^, then 

Vx G Uy^h G ’M : x + hv e U. (4) 

Using this language, the new statement of the Fermat-Reyes theorem is the 
following. The above notation y^h is very convenient whenever thickenings 
(see [12, Lem. 23].) are implicit, as in the following: 

Theorem 2. Let U he a Fermat open set of ’M, and let f : U —^ ‘M 
be a smooth function. Then there exists one and only one smooth map 
r G such that 

V(x,/i)GU: f{x + h)=f{x)+h-r{x,h) in *M. (5) 

Without citing the thickening U, we can state (5) as 

Vx G G *M : f{x + h) = f{x) + h-r{x,h) in *M. (6) 

We ean thus define f'{x) := r(x, 0) G *M for every x £ U. Moreover, if 
f{x) = *a{p,x), Vx G *U C [/ with a G C°°{P x U,M), then 

f'{x)= in ’M. 

Proof. Take (x, h) G U. We first consider the case °h = 0. Since x G U, we 
can write /jv,, = *axipx, -)|v^, where G C°°{Px x 14,K), Vx '■= *Vx C U 
is a Fermat open neighborhood of x and *Px is a Fermat open neighborhood 
of some fixed parameter px G Set Tx := *Vx for the thickening of *Vx. 

Since Tx is a Fermat open subset of (see [12, Lem. 23]), we can find two 
open subsets ^4, R of M such that *{A x R) C Tx, °x £ A and °h = 0 £ B. 
Since U is also a Fermat open set, without loss of generality, we can assume 
that *{A X B) C U. By definition of thickening, we have a + s ■ b £ Vx 
for all a G 4, 6 G R and s £ [0,1]®. This proves correctly property (22) 
of [12]. From this point on, we can continue like in [12], proving that we 
can define a smooth function r : *(j 4 x R) —>■ 'M that satisfies (5) in its 


domain. If °h ^ 0, then h is invertible, and hence we can simply define 
r{x,h) := in a neighborhood ‘{A x B) C U of {x,h) such that 

0 ^ B. Using the sheaf property in *C°°, exactly like in [12], we have the 


conclusion. 


□ 


Here is a simple application. If : M —>■ M is a smooth function, 
then (*gy{x) = *{g'){x). The Fermat-Reyes theorem hence implies that the 
function M defined by 



is smooth. 

Note that the Fermat-Reyes theorem permits to differentiate smooth 
functions only at interior points (in the Fermat topology). 

3. A uniqueness theorem 

Before discussing integration of quasi-standard smooth functions, we 
prove in this section a very useful uniqueness theorem. 

In order to state the theorem, we need another topology on Fermat reals, 
called the w-topology. We summarize the basics of the w-topology; see [14] 
for more details. 

Let dui : x ’M"’ —)• M>o be defined by duj{x,y) = ||°x — °y\\ + 

— yi). Then d^i is a complete metric on ’M”', and the topology 
induced by this metric is called the oj-topology, which is strictly finer than 
the Fermat topology. For any x G ’M” and any s € (0, 1]r, every element 
in the open ball Bs{x-, d^j) is of the form x r for r G MA with ]|r|| < s (see 
[14, Thm. 7, 8, 11]). We will denote by the topological space given by 
a subset A of ’M”" with the w-topology. Note that for U, V open sets of 
Euclidean spaces, not every ordinary smooth function / : U — V induces 
a continuous function */ : *17^^ —>■ For example, if / : M —>■ M is 

the square function f{x) = then (*/)“^(Ri( df; d^;)) =: is not open in 

In fact, one can check that dt 2 G X, but no w-open neighborhood of 
it is contained in X. 

Remark 3. 

(i) The set of all invertible elements in *M is open and dense in both 
and *M, where the latter has the Fermat topology. 

(ii) The ideal is closed in both and *M. More generally, every 
ideal in *M is closed in but only the (unique) maximal (prime) 
ideal 77oo is closed in *M. 

(iii) The sub-topologies on M” via the embeddings M” —> 'M” and MA — 

both coincide with the Euclidean topology. 
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The following lemma is a trivial consequence of the uniqueness of de¬ 
composition of Fermat reals: 

Lemma 4. Assume that x G *M can be written as x = r + 

where G N, r, ai, ai,..., un, oat G M and oi > ... > oat > 1. Then x = 0 

implies r = 0 and a* = 0 for all i = 1,... ,N. 

Proof. Since r = °x = 0, we have the first part of the conclusion. If = 0, 
the second part is trivial. Otherwise, x = ' dtoi is the decom- 

QiT^O 

position of x, and the conclusion then follows from the uniqueness of the 
decomposition. □ 

A writing of the form x = r -|- ' dtai satisfying the conditions of this 

lemma (compared to the decomposition of x, we allow afs to be 0) will be 
called a quasi-decomposition of x. 

Here is the uniqueness theorem: 

Theorem 5. Let A be an arbitrary subset of ’M"', and let f,g:A —>■ 'M 
be smooth functions. If f{x) = g{x) for all x in a dense subset of A^, then 
/(x) = g{x) for all x G A. 

Proof. We may assume that g is the zero function. For any a £ A there 
exist an open neighborhood V of °a in M"', an open set P C RP, a fixed 
parameter p G *P, and an ordinary smooth function a : P x V —M such 
that /(x) = *a{p, x) for all x G n A. 

Letp = (pi,... ,pp) G and a = (ai,... ,an) G *M"' be the components 
of p and a. Let pi = °pi + Pu =■ °Pi + k, I = ■ ,P, and = 

°am + =■ °o-m + hm, m = 1 ,..., n, be the decompositions 

of these components with ki, hm their inhnitesimal parts. For r G M"" such 
that ||r|| < 1 and a -I- r G *V fl A, we can write 

f{a + r) = *a{p,a + r) = *a(°pi + ki,...,°pp + kp,°xi + hi,... ,°Xn + hn), 

where Xm = + for m = 1,..., n. We show below that the coefficients 

in a quasi-decomposition of f{a + r) are ordinary smooth functions of r. 
From the above decomposition, we have ki G Dq.^^ and hm G . Thereby, 
denoting by [u] G N the ceiling of n G M and setting A: := (hi,..., hp), h := 
(hi, ..., hn) and c := max /=i,...,p i\qii], [f^iml), we have that (h, h) G 

Using Taylor’s formula (2) for the smooth function a at the point (°p, °a-|-r) 
and with increments (fc, h), we obtain 

/(a + r) = 'a{p, a + r) = ^ dja{°p, °a + r)- , (7) 

iGJ 


10 



where J is the set of all the multi-indices j € such that |j| < c. Since 

we have that 

P n p / Ni \^^n/ Mm 

{k,hy =n ■ n =n ■ n 

1=1 m=l 1=1 Vi=l / m=l \t=l 

i.e., a finite product of elements in the ideal Dc, the decomposition of {k, hy 
can be written as 

Kj 

{k, hy = 'jje ■ dt^ ( 8 ) 

e=l 

for some constants Kj € N, 'jje G K and Sje € M>i (depending on the 
hxed datum a and /, of course). For simplicity, we can also use ( 8 ) for 
j = 0 € NP”''"' if we set Kj = jje = 1, sje = 0 and dto = 1- Substituting ( 8 ) 
into (7), we get 

Ki 

f{a + r) = *a{p,a + r) = Y'^9ja{°p,°a + r) ■ ^ ■ jje • dt*.^. 

j&Je=l 

To arrive at a quasi-decomposition of f{a + r), it remains to gather all the 
summands having the same term dtg.^. Set 

L := {sje I J G T , e = 1 ,..., iT/} , 

let {rci,... jWh} '■= L be an enumeration of L such that rci > ... > wh, 
and finally, for each £ = 1,..., H, set 

Ceir) := 

Y °CL + ^) ■ ^ Sje = we , j G J , e = 1,..., Kj 



) 3p+m 
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Then 

H 

f{a + r) = Y Ceir) ■ dt^,^ (9) 

£=1 

is a quasi-decomposition of f{a + r), and each Cy—) is an ordinary smooth 
function of r in a neighborhood of 0 G M”". The assumption implies that 
we can find a real sequence (r(u))^gis} 4- 0 such that o-|-r(u) G *V (1 A and 
f(a + r{v)) = 0 for all u G N. Therefore, Yl^=i = 0 for 

all V. Lem. 4 yields that Ci{r{v)) = 0 for all i, and hence ^(O) = 0 by 
the continuity of Ci{—). Setting r = 0 in (9), we get the conclusion that 
/(a) = 0. □ 

This theorem will be used frequently in the following sections. 
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4. The constant function theorem 


Our first problem concerns the existence and uniqueness of primitives of 
smooth functions / : [a, b] —"M, and hence the starting of an integration 
theory. We will solve this problem in next section, and in this section, we 
first prove a useful instrument for solving this problem, namely the constant 
function theorem. 

In the following, we use simplified symbols like [a, b] also to denote the 
corresponding subspace of *M in the category *C°°, i.e., for ([a, b] ~< *M) (see 
[12] for the notation (X -< Y)). We set ’M := ’M U {+oo,— 00 } with the 
usual operations, and °a = a if o € {± 00 }, even if we will never consider 
intervals closed at —00 or + 00 . 

It is not hard to prove that if a, 6 G M, then 

• {(a, 6 )R} C (a, 6 ) 

• {[a, 6 ]R} g [a,b]. 

For example, xt := a — is equal to a in “M, and hence it belongs to the 
interval [a,b], but x ^ *{[a, 6 ]]r} because x does not map M>o into [a,ft]®. 
Analogously, a + dt G {a,b) but °{a+ dt) = a ^ (a,ft)®, so a + dt ^ 
* {(a, fe)R}. This implies that the interval (a, b), with a, b £ *M, is not open 
in the Fermat topology if at least one of the endpoints is finite, and thus we 
are initially forced to consider derivatives only at its interior points, where 

int {* {(a, 6 )r}} = int {* {[a, 6 ]r}} = int[a, 6 ] = int(a, b) 

= {x G ’M I °a < °x < °b}. 

Here, of course, int (A) denotes the interior of A in the Fermat topology, and 
we have used simplified notations like int [a, 5] = int {[a, 5]}. 

As a second step, we will define the notion of right and left derivative at 
a and 6 , respectively, at the end of this section. 

The following lemma is needed for the constant function theorem: 

Lemma 6 . Let J be an interval of *M. Then both the interval J and its 
interior int{J) are connected in the Fermat topology. 

Proof. Let a, b be the endpoints of the interval J. We proceed for a, b 
finite and for the interval (a, 6 ), because the proof is similar for the other 
cases. Assume that there exist two non-empty relatively Fermat open sets 
Ui, U 2 of (a, 6) such that {a,b) = C/i U C /2 and C/i n 1/2 = 0. Set I := {°x | 
X G (a, 6 )}. It is not hard to see that I = [°a,° 6 ]R ^ 0. Moreover, setting 
Vi := {°x I X G Ui} for i = 1,2, we have that Vi is a non-empty relatively 
open set of I, / = lA U V 2 and Vi n V 2 = 0) which implies the conclusion 
since I is connected in M. The proof for the case of the interior is similar 
since int(a, 6 ) = {x G ‘M | °a < °x < °b}. □ 
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Here is the constant function theorem, which is a stepstone for the exis¬ 
tence and uniqueness of primitives for smooth functions. 

Theorem 7. Let J he a non-infinitesimal interval of *M. Let f : J —*M 
he a smooth function such that f'{x) = 0 for each x G int{J). Then f is 
constant on J. 


Proof. Let a, b be the endpoints of the interval J. By the uniqueness Thm. 5, 
it is enough to show that / is a constant function when restricted to int( J) = 
*{(°a,°6)R}, and the assumption that °a < °b (non-infinitesimal interval) 
implies that int(J) 0. Let x G int(J), so that 

V°/i G *M : X -|- h is an interior point of (a, b). (10) 

From the Fermat-Reyes Thm. 2 we get 

V°/i G *M : f{x h) = f{x) -|- h ■ f'[x, h], (11) 


Let us say that (10) and (11) hold for |/i| < r G M>o (here r depends on x), 
and fix such an h. Without loss of generality, we may assume that \ht\ < r 
for each t > 0, where {ht)t is a little-oh polynomial representing h. 

For simplicity, set y := f'\x,h]. As proved in Thm. 2 (or Thm. 22 in [12]), 
we can always find an ordinary smooth function a and a parameter p G 
such that for all s G [0, 1]r 

yt = ^ d 20 i{pt-, xt + s- ht) ds V° t G M>o (12) 

Jo 

f'{x -)- sh) = *{d 20 i){p, X + sh) in *M, (13) 


where f denotes the classical Riemann integral. As a motivation, we can 
say that we would like to consider the function y = *{d 2 Q){p,x -|- sh)ds 

with p,x,h G *M. For fixed p,x,h G *M, we understand it as yt above, 
i.e., using little-oh polynomial representative for each Fermat real, so that 
we get an ordinary integral. Note that this integral is independent of the 
representatives we choose. But \sh\ < r, so that x sh is still an interior 
point, and hence by assumption f'{x + s - h) = 0, i.e. *{d 20 i){p, x + s-h) = 0 
in *M. Written in explicit form, this means 

.. d2a{pt,xt +s ■ ht) ^ 

lim - = 0. 

t^o+ t 


From (12), and using the Lebesgue dominated convergence theorem, we have 


1- yt 

hm — 
i^0+ t 


lim - • / d 2 a{pt, xt + s ■ ht) ds 
i^-o+ t Jo 


r lim "» + »■'■<) d, = 0. 

Jo i->0+ t 
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So y = f'[x,h] = 0 in 'M, and hence (11) yields f{x + h) = f{x) for all 
h such that \h\ < r. This proves that / is locally constant in the Fermat 
topology, and from Lem. 6 the conclusion follows. □ 

We finally introduce the notion of left and right derivatives by means 
of an extension of the Fermat-Reyes theorem from a Fermat open set like 
int(J) to the whole interval J. We first need the following: 

Lemma 8 . Let a, 6 € ‘M and let J be a non-infinitesimal interval of *M 
with endpoints a, b. Let f : J —*M be a smooth function. Then there exist 
S E M>o and a smooth function f : {a — 6, b + 6) —> *M such that f\j = /. 

On the other hand, the smooth function x E {x E *M | 0 < °x < 1} i—>■ 
i E *M cannot be extended to a smooth function defined on any interval 
containing an infinitesimal. 

Proof. By dehnition, we can write the function / as a parameterized exten¬ 
sion of an ordinary smooth function both in a Fermat open neighborhood 
of a -|- dt and of 6 — dt, i.e.: 

f{x) = *a{p,x) \/x€*Ar\J (14) 

f{x) = *fi{q,x) Vx E n J, (15) 

where A, B are open sets in M such that °a ^ A and °6 E R, a E C°°{P x 

A,R), fi E C°°{Q X B,R), p €*P C ‘RP, and q&*Q P ’ML Therefore, 

any two of the smooth functions *a{p, —) : * A —y M, /|int(j) and */3((?, —) : 
*B —>■ M are equal in the intersection of their domains. Since °a a A and 
°b E B, we can always hnd some 5 E M>o such that {*^, int(J), *5} is a 
Fermat open covering of {a — 5,b-\- 5) (let us note explicitly that in this step 
we need to use the trichotomy law of the total order < on "M). The sheaf 
property of the space {a — 5,b -\- 6) yields the existence of a smooth 
function / : (a — J, 6 -|- 5) —S' M such that f\j = f\j. □ 

The notion of thickening is defined only for Fermat open sets U] see 
[ 12 ]. A simple way to extend this notion to arbitrary subsets is the following: 

Definition 9. Let C C and v E *M'^. We say that T C x ’R is a 
thickening of C along v if there exist a Fermat open set U C •R'^ such that 
ecu and T = {(x, h) E Lly | x E C}. 

Here is an extension of the Fermat-Reyes Thm. 2: 

Lemma 10. Let a, 6 E *R and let J be a non-infinitesimal interval of *R 
with endpoints a, b. Then the Fermat-Reyes theorem also holds if we replace 
the Fermat open set U with the interval J in the statement of Thm. 2. More 
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precisely, if x & J and a: ~ a (i.e., °x = °a), so that x G J \ int(J), then 
equation (6) shall he replaced hy 

V°/i € *M>o : f{x + h) = f{x) + h-r{x,h) in *M. 

Analogously, use V°/i € *M<o if x c^b. 

Remark 11. Let r : T —*M be an incremental ratio of a smooth function 
/ : J —‘M defined on a non-infinitesimal interval J of *M, and let s : S' — 
*M be an incremental ratio of the restriction /|int(j)- Then s is essentially a 
restriction of r in the sense that 

Vx G int(J) V°/i € *E : r{x,h) = s{x,h). 

We can hence define f'{x) := r{x, 0) for all x G J, obtaining an extension of 
/' from int( J) to the whole J. 

Proof. Let / : {a — 6, b + 6) —*M be an extension of / : J —*M as 
in Lem. 8, set U := int(a — (5, 6 -|- <5) and let f'[—,—] ■ U —)• *M be the 
incremental ratio of /. Then T := |(x, h) G [/ | x G j| is a thickening of 

J and r := f'[—, — ]|t is a searched incremental ratio of /. The uniqueness 
part follows from the uniqueness Thm. 5. □ 

5. Existence and uniqueness of primitives 

We can now prove the existence and uniqueness of primitives for smooth 
functions defined on an interval [a, b] with finite endpoints. 

Theorem 12. Let a, b G with °a < °b, let f : [a, b] —*M be a smooth 
function, and let u G [a, b]. Then there exists one and only one smooth 
function 

I : [a, b] —^ ‘M 

such that 

l'{x) = /(x) Vx G [a,h], and 
I{u) = 0. 

Proof. We may assume that u = a, since if I' = / on [a,b] and I(a) = 0, 
then J(x) := I{x) — I{u) verifies J' = I' = f on [a, b] and J{u) = 0. 

For every x G [a,b], we can write 

f\Vx— ~)\Vx 

for suitable px G *Ux C •RP'', Ux an open subset of 14 an open subset 
of R such that x G *14 H [a,b] =: Vx and Ux G C°°(Ux x 14,1^)- We may 
assume that the open sets 14 are of the form 14 = (°x — dx,°x -|- for 
some 6x G R>o- 
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The idea for the construction of the function I is to patch together suit¬ 
able integrals of the functions *ax{px, —) and, at the same time, to respect 
the condition I (a) = 0. In order to get a quasi-standard smooth function, 
we have to patch together integrals in an order so that each one is the “con¬ 
tinuation” of the previous one. In other words, on the non-empty connected 
intersection of the domains of any two of these integrals, the integrals must 
have the same value at one point, so that we can prove that they are equal 
on the whole intersection. 

Since is an open cover of the compact set [°a,°6]]R, we can 

hnd a finite subcover. That is, we can find xi,... ,Xn G [°a,°6 ]]r such that 
(I4i)j=i n open cover of [°o, °6 ]r. We will use simplified notations 

like .— ^Xi: etc. 

By shrinking each Vi if necessary, we may always assume to have chosen 
the indices z = 1 ,..., n and the radii Si € M>o such that 


°a = xi < X 2 < ■ ■ ■ < Xn = °b, and 


Si Xi-^\ Si-^\ ^ Xi Si “t” Vz — 1, . . . , U I. 

In this way, the intervals Vi = {xi—Si,Xi+Si)s. and V^+i = (xi+i—Jj+i, Xj+i-|- 
intersect in (xj+i — dj+i, Xi-|-5i)R. For each z = 1,... , zz — 1, we choose 
a point Xi G (xj+i — Si+i,Xi -|- (5i)R and define, recursively: 



(16) 


Jy 

/i(x) : = *zi(pi,X,a) Vx G *Tl. 


X 



ii+i{q,x): = f ai+i{q,s)ds Vg' G 14^ Vx G T^+i Vz = 1,..., n - I, (17) 


4+i(x) : = •ij+i(pi+i,x) -I- Ii{xi) Vx G 'Vi+iMi = 1,... ,7Z - 1. 

So every 4 is a quasi-standard smooth function defined on* Vi, and moreover, 
from Thm. 2 it follows that 

I'i{x) = *ai{pi,x) = f{x) \/xe*Vi. 

Therefore, for each point x in 

•14 n ’ 17+1 = 'iv n Vi+i) = ’{(xi+i - Si+i,xi + 4)m} 

= int([xi+i - Si+i,Xi -k 4]), 


we have 


I'iix) = fix) = I'i+iix), and 
(4 - 4+i) (xi) = 0. 


(18) 
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So, from the constant function Thm. 7 it follows that = /j+i on [xj+i — 
5i+i,Xi + (5i]. We can hence use the Fermat open cover {*Vi fl [a,b])-^^ ^ 

of the space [a, b] G to patch together the functions to obtain 

a smooth function I : [o, b] —*M, by the sheaf property of smooth func¬ 
tions. This function satishes the required conditions because of (18) and 
the equalities I {a) = Ii{a) = 0. 

To prove the uniqueness, suppose that J verifies J' = / on [a, b] and 
J{u) = 0. Then again by the constant function Thm. 7, we know that 
(J — I)\[a,b] is constant and equal to zero at u, and hence a zero function. □ 

The above theorem can be extended to smooth functions defined on an 
arbitrary non-inhnitesimal interval. 

Theorem 13. Let a, 6 € *M and let J be a non-infinitesimal interval o/*M 
with endpoints a, b. Let f : J —> ’M be a smooth function, and let u & J. 
Then there exists one and only one smooth function L : J —)■ ’M such that 
I' = f and L{u) = 0. 

Proof. We proceed for the interval (—oo,+oo), since the other cases are 
similar. Set r := °u, and for each k € N>o, dehne 


fk • f\[r—k,r+k]- 

By Thm. 12, for each k there exists a unique smooth function : [r — fc, r + 
k] — > *M such that I^ix) = fk{x) = f{x) for x € int[r — k,r + k] =: 14 
and Ik{u) = 0. Note that (14)fc>o i® a- Fermat open cover of ’M. Moreover, 
since I'f.{x) = f{x) = I'j{x) for every x G 14 n V)- and Lk{u) = Ij{u), by 
the constant function Thm. 7, 4 and Ij coincide on 14 H V)'. By the sheaf 
property of smooth functions, we get 

3!/g’C°°CM,’K) : I\v,=Ik VfcGN>o. (19) 

Note that for any x G *M, there exists k = k{x) G N>o such that x G int(14). 
Therefore, for each x G *M and for /i G *M sufficiently small, we have 

L{x + h) = 4(x + h) = L{x) + h ■ L'[x, h] = 4(x) + h ■ L'fiyx, h]. 

This and the uniqueness Thm. 5 imply that the incremental ratios of I and 
4 are equal, i.e., L'[x,h] = L'j^[x,h]. Thus, I'{x) = I'j^{x) = /(x), and finally 
I{u) = Ii{u) = 0. 

This proves the existence part. 

Since the restriction of a primitive of / to [r — k,r k] is a primitive of fk 
for each k, the uniqueness then follows from Thm. 12. □ 

We can now dehne 
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Definition 14. Let J be a non-infinitesimal interval of *M. Let / : J —>• *M 
be a smooth function, and let u & J. We define ^ / to be the unique 
function verifying the following properties: 

r(-) r(-) 

(i) / f '— /(®) ^ is ^ smooth function; 

J U J U 

pu 

(ii) / / = 0; 

J U 

(hi) Vx € J : {^j {x) = ^ J f{s) ds = f{x). 

Remark 15. 


1. Note that, although we require the interval J to be non-infinitesimal 
for the domain of the smooth function /, the integral / is also 
defined when u € J and u ~ u. 

2 . It is not hard to see that if /, 5 : J —)• *M are smooth functions, 
a,be J with a<b, and f\^a,b) = 9\{a,b), then f = g. We get the 
same conclusion if °a < °b and we assume only /|int(a,fe) = 5lint(a,fe)- 

3. From Def. 14, we obtain a generalization of the usual notion of integral. 
Indeed, for a, 6 , u € M with a < u < 6 , let / € C°°([a, 6 ]r, M). We 
extend / to an ordinary smooth function defined on an open interval 
{a — 5, b + (5)r, for some 6 G M>o. We still use the symbol / to denote 
this extension. We consider the quasi-standard smooth function 

f{s) ds^ :*{a-6,b + (5)r — > 'M. 

Now, we have that 

[a,b] C ’(a - (5 )r 

so that we can consider the restriction I\[a,b]- It is not hard to prove 
that this restriction verifies all the properties of the previous Def. 14 for 
the function */. Meanwhile, because I is the extension of an ordinary 
smooth function, it also verifies 

f{s) ds G M. 

( 20 ) 

The recursive definitions (16) and (17) applied to f\[u,v] yield the following 
useful result: 


Vx G 


PX ® / PX \ 

J *f{s)ds = I{x)= N- /(s)dsj = 
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Lemma 16. Let a, 6 € *M and let J be a non-infinitesimal interval of *M 
with endpoints a, b. Let f : J —’M be a smooth function, and let u, v ^ J 
with u < V. Then there exist n G N>o, y, y, 6 G M”, {pi)^=i such 

that 

(i) Ui G C°°{Ui X yi,E) and pi G *Ui, where Ui C RP^, Vi = {yi - 6i,yi + 
(5i)K, 5i G M>o, \/i = l,...,n 

(a) f{x) = *ai{pi,x) Vx G *Vi n J Vi = 1,..., n 

(Hi) yi G {yi+i - 5i+i, y* + (5i)R = Vit^ V+i / 0 Vi = 1,..., n - 1 

(iv) °u = yi < y 2 < ■ ■ ■ < yn = °v if °u < °v and °u = xi = °v otherwise 

(^) In f = /r *ai{pi,s) ds + * 02 (^ 2 , s) ds + ... + ’an{Pn, s) ds, 

and a little-oh polynomial representing f is given by 

ryi rt>t 

tGR>o>-^T ai{{pi)t,s)ds+ T Q;n((pn)t, s) ds G M. 

Jut Jyn-l 


6. Classical formulas of integral calculus 

In this section and the first part of next, we will prove the classical 
formulas of integral calculus. 

The first results concerning the integral calculus are the following: 

Theorem 17. Let J be a non-infinitesimal interval o/*M, let f,g:J —>■ *M 
be smooth functions, and let u, v G J. Then we have: 

(^) /:(/+v)=/:/+/:v 

/:a/=a/:/ vag-m 
(Hi) fu f = f(.v) - f{u) 

luif ■ a) = [f ■ aVu- fuif ■ a') 

We will prove the additivity with respect to the integration intervals, i.e., 
lu ■f /iT f ~ lu after the proof of the commutativity of differentiation 
and integration in next section. The latter relies on the quasi-standard 
smoothness of the function x jQ^/(s,x)ds, which will be proved as a 
consequence of the *C°° smoothness of suitable infinite dimensional integral 
operators. 

Remark 18. In fact, one can prove the above properties together with the 
additivity with respect to the integration intervals using Lem. 16. As an 
example, we prove in this remark that /J / + / for any smooth 

function / defined on a non-infinitesimal inteval J of *M with u,v,w G J. 
By Lem. 16, we can write 
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I-V ryi 

/ /=/ *aiipi,s)ds + . 

rv 

■■+ *aniPn,s)ds 

J u J u 

Vn—l 

nw rzi 

nw 

/ /=/ * fdiiqi, s)ds + .. 

J V J V 

Zm — l 


Since the function / is quasi-standard smooth, we may assume that «„ = /3i, 
i.e., an{pi-,x) = f{x) = /3i(gi,x) for every x € (y„_i,zi). The conclusion 
then follows from the additivity of ordinary integral with respect to ordinary 
integration intervals. 

This proof uses local expressions of quasi-standard smooth functions. 
The proofs given below are of “global” nature, which is the point we try to 
emphasize in the present work. 

Proof. In the following, we will use the notations 



To prove (i), it suffices to note that F + G is smooth and 
(F + G) (n) = F{u) + G{u) = 0, 

{F + G)' (x) = F'{x) + G'{x) = fix) + gix) = [f + g) (x) Vx G J. 

From the uniqueness of Def. 14, the conclusion F + G = ^ (/ + 5 ) follows. 

For the proofs of the other properties, we can follow the same method: 
we define a suitable quasi-standard smooth function H starting from the 
right hand side of the equality we need to prove, and we prove that H 
verifies the same properties that characterize uniquely the function on the 
left hand side. For example, to prove (ii) we consider i/(x) := AF(x); to 
prove (iii) we define H(x) := fix) — fiu): and finally to prove (iv) we set 
Hix) := fix) ■ gix) - fiu) ■ giu) - f^if ■ g'). □ 

Here is the result of integration by substitution: 

Theorem 19. Let J, Ji be non-infinitesimal intervals o/*M, let u, v £ Ji, 
and let 

Ji J ’E 


be smooth functions. Then 




p'is) ds. 
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Proof. Define 


F{x) 

Hiy) 

G{y) 


f f VxeJ 

J <~p{u) 

f‘p{y) 

/ / Vy G Ji 

J ip{u) 


Vy G Ji. 


Each one of these functions is smooth because of Def. 14 (we recall that u G 
Ji is fixed) or because it can be written as composition of smooth functions. 
We have H{u) = G{u) = 0, H{y) = F [(/?(y)] for every y G Ji, and by the 
chain rule ([12, Thm. 29]) H'{y) = F’[ip{y)] ■ (p'{y) = f [(p{y)] ■ (p'{y) = G'{y) 
for each y G Ji. From the uniqueness of Thm. 12 the conclusion H = G 
follows. □ 


7. Infinite dimensional integral operators 

If we want to prove, in our framework, that 

^(1 ( 21 ) 

we need to prove that the function defined by 

xG[a, 6 ]i-7> f /(s,x)dsG*M (22) 

Jo 

is smooth. Of course, in (21) |^(s,x) := [/(s, —)]^ (x). The smoothness 
of (22) follows as a trivial consequence of the smoothness of the function 
defined by 

4' : if,u,v) G •C°°([a,6],*M) x [a,b]^ ^ T / G ‘M. (23) 

J U 

In fact, the function defined in (22) can be written as /g(—,0,1) o (/ o r)^, 
where^ 

r : (x, s) G [a, b] x [0,1] )->■ (s, x) G [0,1] x [o, b]. 

Therefore, the smoothness of (22) follows from the smoothness of (23) and 
the Cartesian closedness of the category of Fermat spaces. For simplic¬ 
ity, in the following we will denote by the space *C°°([o,6],‘M). 


®See [12, Def. 1] for the notations (—)^ and (— 
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From Thm. 19, we have 


rv pi 

/ / = / f [y ■ - u) + u] ■ {v - u) dy. 

J u Jo 

Thereby, setting 

V’(/, u,v) :y e [0, l]i-^ f[y-{v-u) + u]-{v-u) 

I ■. f e ^ [ / G ’M, 

Jo 

it is easy to show that 

'll; : X [a,bf is smooth 

and I^{f,u,v) = I [ip{f,u,v)]. Therefore, the following result suffices to 
prove that the function is smooth. 

Lemma 20. The function 

life ^ / / G ’M 

Jo 

is smooth, i.e., it is an arrow of the category *C°°. 

Proof. Take a hgure of i.e., let 

S C ’R" 

dG •C°°(5 ,*r[°’^]). (24) 

Recall that here we use the simplihed symbol [0,1] to denote the correspond¬ 
ing subspace of ’R, i.e., ([0,1] ^ ‘R). From Property 5(i) of Thm. 18 in [12] 
we have that [0,1] = ([0,1] ^ *R) = [0,1]. Therefore, from (24), we know 
that dP : S X [0,1] —?> *R is smooth. We need to prove that the function 

s G 5 !->■ / d{s) = f d^(s, u) du G *R (25) 

Jo Jo 

is smooth, i.e., that it belongs to S*R°°(S', *R). From Properties 5(h) and 
5(i) of Thm. 18 in [12], we have that 

S X ((blj ={S < ’R") X ([0,1] ^ ’R) 

= (^{S X [0,1]) ^ 

= Sx [0,1]. 

Hence, dy : S x [0,1] —*R is smooth, i.e., 

dy : S X [0,1] —> *R is an arrow of S*R°°, (26) 
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because S*M°° is fully embedded in the category of Fermat spaces (see 
[12, Thm. 18(4)]). Property (26), for every {s,u) € S' x [0,1], yields the 
existence of 

Us,u Fermat open neighbourhood of ps^u defined by 

Vs,u Fermat open neighbourhood of (s,tt) dehned by Vs^u in S x [0,1] 

as,u G C°°{Us,u X P.,n,K) 

such that d^\vs,u — *ces,u{Ps,u, —)\Vs,u- Now we fix s G S, and everything 
above only depends on u € [0,1]. Since (s,u) € = *14 n (S x [0,1]) and 
14 is open in we have that (°s, °u) € 14- Thus, we may always assume 
that 14 = Au X Bu, for some open subsets Au in and in M, with 
(°s,°tt) e Au X Bu- Hence, (s,tt) G * {Au x Bu) = *Au x *Bu = *14- In 
this way, we obtain an open cover {Bu)u£[o i]i; [!*’ I]®’ from it we can 
extract a hnite subcover. Therefore, we can hnd ui,... ,Un G [0, 1]r such 
that {Bui)^^i open cover of [0, IJr. As in Thm. 12, we may assume to 
have chosen the indices i = 1, ..., n and have found suitable 4’s in M>o so 
that 


0 = Ui < U2 < ■ ■ ■ < Un = I 
ddui — {Ui di,Ui (5j)R 

ddui fl — {ui+i 4+1) Ui + 4)r 7^0 — l,...,?r 1. 

Finally, set A := Au^ fl • • • fl Au„, which is an open neighbourhood of °s, so 
that s G *A. Take a point Ui G (tti+i — 4+i) 'Wi + 4 )k- For every x G *A fl 5, 
we have 

dA (x,u) = ’aui {pui ,x,u) Vn G n [0,1], 
and hence Lem. 16 implies that 

rui ru2 

/ d^{x,u)du = / *aui{pui,x,s)ds + / *au 2 {Pu 2 ,x, s) ds + ■ ■ ■ 

Jo Jo J ui 

+ / *auu{Puu,x,s)ds =:*'y{p,x), 

J Un—1 

where p := (Pui, ■ ■ ■ ,Pun) G Um x ••• x Uu„ 4 . From the 

classical version of the differentiation under the integral sign, we have that 
7 G C°°{U X A,M), where U := C4i x ••• x C4„, and hence we get the 
conclusion. □ 


Corollary 21. Let J be a non-infinitesimal interval of *M . Then the func¬ 
tion 

PV 

if, u, v) G •C°°( J, ’M) X J2 ^ / / G ’M 

J U 

is smooth. 
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We can now prove the differentiation under the integral sign: 


Lemma 22. Let J be a non-infinitesimal interval o/*M, and let f : [0,1] x 
J —> *M be a smooth function. Then 

Proof. From Cartesian closedness, it follows that the map x G J /(~5 x) G 
•]^[0d] jg gixiooth. Together with Lem. 20, we know that 

X G J [ f{s, x) ds G *M 

Jo 


is smooth. Let d 2 f [—be the incremental ratio of / with respect to its 
second variable. Then for each x G int(J) and for /i G *M sufficiently small, 
we have 


[ f{s,x-\-h)ds= [ {f{s,x)-\-h-d 2 f[s,x;h]}ds 
Jo Jo 

= f f{s,x)ds-\-h- f d 2 f[s,x-,h] 

Jo Jo 


ds. 


Therefore, the uniqueness of the incremental ratio R[—, —] of the function 
fo /('S) ~) bs yields that 


Hence, 


R[x,h]= f d 2 f[s,x;h]ds V°/i. 

Jo 


4- ( [ f{s,x)ds ) = i?[x,0] 


dx 


J d 2 f[s,x; 0 ]d 


as expected. 


□ 


The following lemma is needed for the “global” proof of the additivity 
of integral with respect to the integration intervals; 

Lemma 23. Let J be a non-infinitesimal interval of *M, let f : J —*M be 
a smooth function, and let v G J. Then 


dx 



= -fix)- 
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Proof. Set (p{s) := x + {v — x)s for each s G [0,1]. From Thm. 19 (integration 
by substitution), we get 

rv rl 

/ f{t) dt= f [x + {v — x)s] • (u — x) ds. 

Jx Jo 

From this, using Lem. 22 (differentiation under the integral sign), the prod¬ 
uct rule ([12, Thm. 28]) and the chain rule ([12, Thm. 29]), we get 


_d_ 

dx 


f^ = - f[x + {v-x)s\ds+ 


+ {v - x) 


— ~ f [x + {v — x)s] ds+ 

Jo 

+ [ [ ,s-f'[(p{,s)]-(p'{,s)d,s. (27) 

Jx Jo 

Integration by parts (Thm. 17(iv)) applied to the last integral yields 

[ s- f [v7(s)] • ip'{s) ds = [(/ o (^)(s) • s]o - / (/ o </5)(s) ds 
Jo Jo 

= f{v) - [ f[x + {v- x)s] ds. 

Jo 

Substituting in (27) and using Thm. 17(iii), we obtain 




[x -b (u — x)s] • (1 — s) ds 


dx 


+ / f [x + {v — x)s] ds 

Jo 

= -fix), 


as desired. 


□ 


Now we are led to the “global” proof of the additivity of integral with 
respect to the integration intervals: 


Corollary 24. Let J he a non-infinitesimal interval of *M, let f : J —>• *M 
he a smooth function, and let u, v, w ^ J. Then we have: 


rv nw rli 

J U J V J u 


/■ 


Proof. Set F{x) := fufJ-f^^f for all x € J. Then Lem. 23 yields that 
F'(x) = /(x) — /(x) = 0 at each x G J. Hence by Lem. 6, F is constant and 
F(v) = F(u}, which is our conclusion since f = 0- Cl 
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Using Cor. 24, we can calculate the incremental ratio of the integral 
function ^ f and get a form of mean value theorem called Hadamard’s 
lemma as follows: 

Corollary 25. Let J he a non-infinitesimal interval o/*M, let f : J —>■ *M 
be a smooth function, and let u £ J. Then 

(i) for each x £ J and h £ ‘M, if x -\- h £ J, then f{s)ds = 
h ' Jq f{x + hs) ds. Therefore, the incremental ratio of the smooth 
function ^ f is 

[x,h]= J f{x-\-hs)ds Vx G int(J) V°/i € *M; 

(a) Vx G int(J) V°/i G “M : f{x-\-h)—f{x) = h-f^f'{x-\-hs)ds. 

Proof, (i): The first statement follows from Thm. 19, and the second state¬ 
ment then follows from Cor. 24 and Thm. 2. 

(ii) : It follows by applying (i) to the smooth function f' : int(J) —*M 

together with Thm. 17(iii). □ 

Standard and infinitesimal parts of an integral 

To find a useful formula for the standard and the infinitesimal parts of 
an integral, we introduce the following: 

Definition 26. Let S C and let / : 5 —be a smooth function. 
Then 

(i) °S := {°x G I X G 5}. 

(ii) If °S C S, then we set °/ : r G °5 i->- ° [f{r)] £ which is called the 
standard part of /. 

(hi) Assume that °S is an open subset of M", and * {°S) = S. Then we set 
V •= * (°/) and 6f := f — ^f, which are called, respectively, the shadow 
and the infinitesimal part of the smooth function /. Analogously, we 
set 5x := X — °x for any x G *M. Note that '^f{r) = °/(r) for each 
r G °S. 

For example, if S := int(a, 5), then °S = (°a, °6 )r C S, °S is an open 
subset of M, and * (°5) = S. Note that the standard part °/ is a function 
defined and with values in standard points, whereas the shadow ’^f can 
have nonstandard values for nonstandard points in the domain. Therefore, 
{6f){x) is not necessarily equal to 6{f{x)), unless x £ °S C S] see (i) of the 
following lemma. The correctness of the dehnition of shadow is proved in 

(iii) of the following lemma: 
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Lemma 27. Let S C ’M”. Assume that °S C S, and let f : S —be a 
smooth function. Then 

(i) =°[f{x)] VxeS'. 

(a) If°S is an open subset ofW^, then °f € 

(Hi) If°S is an open subset o/M"’ and * {°S) = S, then “/ € *C°°(5,’M'^) 
and 6f € L>^). 


(iv) Under the assumptions of (Hi), if we further assume that n = 1, 
S D int(a, b) and u, v £ int(a, b), then the standard part of the integral 
fu f given by 


O 



7 



and the infinitesimal part is given by 


nv \ n°u rv nv 

f)= 7 + / 7 + / 7 

J u / J u J°v J u 


V 

= E ^ f h svY-Y. ‘ f h suy+ I SI 

i=l *• i=l 

S Doo, 


j! 


f 


where [a;(tt)] denotes the integer part ofu}{u). 

As trivial consequences of (iv), any integral over an infinitesimal inte¬ 
gration interval is infinitesimal, and any integral of an infinitesimal-valued 
smooth function is infinitesimal. 

Proof, (i): We can write f{y) = *a(p,y) for each y € *V D S, where 
°x £ V and V is an open subset of M”. Therefore, ( 7 ) (°a^) = ° [f{°x)] = 
° [*a{p, °x)] = a{°p, °x). On the other side, ° [/(x)] = ° [*a{p, x)] = a{°p, °x). 

(ii) : We write again = *a{p, where now V is an open neigh¬ 

borhood of r £ °S C S. But y n °5 C *V O S, so 7lyn°S ~ “(7; “)l\/n°S’ 
y n°5 is an open neighborhood of r and a{°p, ^ smooth function. 

This proves our claim. 

(iii) : This is clear from the property of the extension *(—) of ordinary smooth 
functions and from the assumption * {°S) = S. Observe that, using the 
previous notations, ^f{x) = * (°/)(x) = *a{°p,x) for x £ S (with a and p 
depending on x). Hence, {6f){x) = /(x) — “/(x) = *a{p,x) — *a{°p,x) £ 

Dd 
^oo • 

(iv) ; Lem. 16 yields 


O 



O 



*Q!i(pi, s) ds + ... + 



‘an{Pn,s)ds 
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The inside of the brackets of the right hand side of this equality is a smooth 
function of u, xi,..., Xn-i, v, pi,... ,pn- Taking the standard parts (that, 
as usual, includes the use of the dominated convergence theorem for f —>■ 0 ^ 
in M) and using ( 20 ), we obtain that 


( nv \ pxi n°v 

/ /) = / ai{°pi,s)ds + ... + ani°Pn,s)ds 
Ju J J°U Xn —1 


Xi 


ai{°pi,s)ds +... + + an{°Pn,s)ds 
Jxn-1 


7- 


Now, we can compute the infinitesimal part by 


/ = 


PV P°V 

f- 7 

J u J °u 

P°U P°V PV P°V 

/ (7 + 7)+/ (7 + 7)+/ (7 + 7)-/ 7 

Ju J°u J°v J°u 

p°u p°v-\-6v PV 

I 7+/ 7+/ 7 

J °v Ju 

{dvT + r 

J u 


' °u-\-5u J ^v 




7, 


where, in the last step, we have used the infinitesimal Taylor’s formula (2) 
for (the extension of) ordinary smooth functions. □ 


The standard and infinitesimal parts of general Fermat spaces will be 
discussed in [15, Sec. 4.2]. 


Example: Divergence and curl 

Now we use suitable infinitesimals in Fermat reals to revisit the classical 
concepts of divergence and curl. 

Traditionally, divj4(x) is understood as the density of the flux of a vector 
field A € C°°(C/,M^) through an “infinitesimal parallelepiped” centered at 
the point x in an open set U C M^. To formalize this concept, we take three 
vectors hi, h 2 , 7 € and express them with respect to a fixed base ei, 
62 , 63 G as 


hi = k} ■ 61 + kf ■ 62 + fef • 63 , where kj € ’M. 
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We say that P := 
if 


(x, hi, ^ 2 , / 13 ) is a (third order) infinitesimal parallelepiped 
X € and each kj E D 3 . 


Here the requirement that each E D 3 is to make sure that the multipli¬ 
cation of any four (and hence more) such kj ’s is zero. The flux of the vector 
field A through such a parallelepiped (toward the outer) is, by definition, 
the sum of the fluxes through every “face”: 




X — -/13 + thi + sh2) • (^2 X hi) ds+ 


+ j dt [ A{x + ^h 3 -I- thi + sh2) • (hi x (12) ds -|- 


where the fluxes through the two opposite faces spanned by hi and h 2 are 
explicit in the above formula, x is the cross product, and the dots • • • 
indicate similar terms for the other faces of the parallelepiped. Let us note 
that e.g. the function s A{x — ^h 3 +thi + sh 2 ) is a quasi-standard smooth 
function. In this case, the hxed parameter is p = (x, t, hi, h 2 , hfi) E It 

is easy to prove that if H E C°°{U, M^) and the oriented volume Vol(hi ,h 2 ,hfij 
of the inhnitesimal parallelepiped P = (x, hi, h 2 , hs), is not zero, then the 
following ratio between infinitesimals exists and is independent of hi, h 2 , 
(see [9] for the notion of ratio between infinitesimals): 

divH(x) :=-^— • [ A • n dS gR. 

Vol(hi, h 2 , h 3 ) Jp 

To define the curl of a vector field A E we say that C := 

(x,hi,h 2 ) is a (second order) infinitesimal cycle if 

X G U and each kj G 02- 

Here the requirement that each kf G D 2 is to make sure that the multiplica¬ 
tion of any three (and hence more) such kPs is zero. The circulation of the 
vector field A on this cycle C is dehned as the sum of the “line integrals” 
over each “side”: 

I A-t dl := I A{x — h 2 +thi)-hi dt— j A{x+-h 2 +thi)‘hi dt+-■ ■ , 

Jc J-\ 2 J_i 2 

where the line integrals of the two opposite sides spanned by hi are explicit 
in the above formula, and the dots • • • indicate similar terms for the other 
sides of the cycle C (noting the consistency of the orientation on each side). 
Once again, using exactly the calculations frequently done in elementary 
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courses of physics, one can prove that there exists one and only one vector, 
curl74(x) € M^, such that 

/ A • t dl = curM(x) • (hi x / 12 ) 

Jc 

for every infinitesimal cycle C = (x,/ii,/i 2 ), representing thus the (vector) 
density of the circulation of A. 


8. Inequalities for integrals 

In this section, we derive similar inequalities for integrals of smooth func¬ 
tions to the classical ones for integrals of ordinary smooth functions. Thanks 
to the total order on *M, we are able to prove the monotonicity of integrals, 
a replacement of the inequality for absolute value of smooth functions, and 
the Cauchy-Schwarz inequality, with respect to arbitrary integration inter¬ 
val. The strategy is, we hrst prove inequality for real integration interval, 
and then use integration by substitution for arbitrary integration interval. 

Recall (see e.g. [10]) that if x, y € *M, then x < y if and only if for any 
xtiVt £ Ko[i] with X = [xi] and y = [yt], there exists zt G Mo[t] such that 
[zt] = 0 in ’M and x* < -|- zt for t G M>o sufficiently small. In proving 

these integral inequalities, we will use [10, Thm. 28]: if {xt)t is any little-oh 
polynomial representing x G *M, then x > 0 is equivalent to x^ > 0 for t 
sufficiently small, or x = 0 in ’E (i.e., xt = o{t) as t —)• O"*" in E). 

The first two results concerns the monotonicity of integrals: 

Lemma 28. Let J be a non-infinitesimal interval o/*E, and let f, g ■ J — 
*E be smooth functions. Let a, 6 G J PI E, with a < b. Then the assumption 

g{x) < fix) Vx G [a,6]R (28) 

implies g < j(f. 

Proof. It suffices to prove the conclusion for y = 0. By Lem. 16, we can 
write 

rb ^ rXi 

/ f = '^ ‘aiipi,s)ds, 

J a J Xi-i 

where ai G C°°iUi x 14,E), Ui C EP% pi G *Ui, I 4 = (x* - 6i,Xi + (5i)R, 
6i G E>o, Xi G (xi+i - 5i+i,Xi -K (5i)R = 14 n 14+1 C (xi,Xi+i)R, and 

fix) = *aiipi,x) Vx G *14 n J, (29) 

a = Xq = Xi < X2 < . . . < Xn = Xn = b- 

For each i = 1,... ,n and each s G [xj_i,Xj]R, we have a = xq < Xj_i < 
s < Xi < Xn = b, so s & [a, 6 ]r. Moreover, Xj_i G 14 and Xj G 14 yield that 
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Xi — 5i < Xi-i < s < Xi < Xi + 5i, so s € Vi r\ J. Thereby, assumption (28) 
and equality (29) imply that f{s) = *ai{pi,s) > 0. The mentioned result 
[10, Thm. 28] yields that 

> 0 ; ai{{pi)t,s) > 0 

or (30) 

‘ai{pi,s) = 0 in *R. 


Set Wi(t,s) := 0 if ai{{pi)t,s) > 0 and Wi{t,s) := ai{{pi)t,s) otherwise. 
Then, 

Vt G M>o Vs € [xi_i,Xj]K : ai{{pi)t, s) > Wi(t, s), (31) 

and from (30) we get that Wi{t,s) = o{t) as t —> 0+ in M. Note also that 
Wiit, —) : [xi_i,Xj]]R —> M is integrable by definition, so that (31) implies 

T ai{{Pi)t,s)ds > j- Wi{t,s)ds Vt G M>o. (32) 

^i—1 ^i—1 


But Lebesgue dominated convergence theorem and the fact that Wi{t,s) = 
o{t) yield 


lim Wi{t, s) ds = rf ' lim ds = 0, (33) 

t —^ 0 + t t —^ 0 + t 


i.e., the right hand side of the inequality (32) viewed as a function of t is 
o{t) as t —)• 0"*“ in M, and hence it is in Mo[i] and it represents 0 in *M. In 
other words, we have proved that 



•ai(pj,s)ds > 0. 


We thus have 


^ rxi 

/ ^ = 

J a i=\ 1 


*ai{pi,s)ds > 0 , 


which proves our conclusion. We finally note that in (33), all the endpoints 
of the integration intervals, which include a = xq and 6 = do not depend 
on t. □ 


We can now extend this result to the case where the endpoints a, b of 
the integration interval are arbitrary Fermat reals. 

Theorem 29. Let J be a non-infinitesimal interval of *M, and let f, g : 
J —)■ ‘M be smooth functions. Let a, b ^ J, with a < b. Then the assump¬ 
tion 

g{x) < f{x) VxG[a,6 ] (34) 

implies f^g <f^f. 
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Note that in this theorem (also the next two theorems), we include the 
case when the integration interval is infinitesimal, i.e., °a = °h. 


Proof. As above, without loss of generality, we may assume g = 0. Using 
integration by substitution (Thm. 19), we can write 

f f = {b-a)-f f[a +(b-a)a]da (35) 

J a Jo 

For each a € [0, 1]]r, we have a < a + {b — a)a < b. Note that a + {b — a)a 
is in general not standard if a is not standard, which justifies the stronger 
assumption (34) compared to the previous (28). Assumption (34) yields 
that / [a + (6 — a)cj] > 0. Lem. 28 implies that the right hand side of the 
integral in (35) is non-negative, from which the conclusion follows. □ 


Although we haven’t justified the validity of integrals of absolute val¬ 
ues of smooth functions, the following result is a possible substitute of the 
property / < /^^ |/| for classical integrals. 


Theorem 30. Let J be a non-infinitesimal interval of *M, and let f, g : 
J —)• 'M he smooth functions. Let a, h & J, with a < b. Then the assump¬ 
tion 


\f{x)\<g{x) Vx€[o, 6 ] 


(36) 


implies 


laf <fa9- 


Proof. The proof is conceptually identical to those of Lem. 28 and Thm. 29. 
We first prove the conclusion for a, 6 € M. The interval [a, 6 ]r can be covered 
by a finite number of intervals W := (x* — 6i,Xi -|- (5i)R, where we can write 
both f{x) = *ai{pi,x) and g{x) = *fii{qi,x) for every x G n [a,6]. 
Furthermore, we have 


pb 

/ / = X] / 'ai{pi,s)ds, 

J CL i = ^ J ^i — 1 

[ 9 = '^f '/3iiqi,s)ds, 

J d ^^ J X-1 _ 


where Xi € (xj+i —(5j+i,Xj-t-(5i)K and a = xq = xi < X 2 < ■ ■ ■ < Xn = Xn = b. 
For each i = 1,..., n and each s G [xj-i, assumption (36) and [10, Thm. 
28] yield that 


G M>o : \ai{{pi)t, s)| < fii{{qi)t, s) 
or 

I'afipi, s)| = *fii{qi, s) in ’M. 
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Set 


0 , a \ai{{pi)t,s)\ </3i{{qi)t,s) 

\ai{iPi)t,s) \ - f3i{{qi)t,s), otherwise. 


Therefore, 


Vt € M>oVs G [a;i_i,Xi]R : \ai{{pi)t, s)\ < I3i{{qi)t, s) + Wi{t, s). (37) 

As in the proof of Lem. 28, one can show that for any fixed s € 

Wi(t, s) = o{t), and for any fixed t G M>o, the function Wi{t, —) : [xj_i, Xj] —>■ 
M is integrable. Again, one gets ^Wi{t,s)ds = o{t), and as a conse¬ 
quence, ^t(;i(t,s)ds G Mo[t] and it represents 0 in ‘M. From inequality 
(37) we get the following inequality of little-oh polynomials: 


rxi 

i=l 


ai{{pi)t,s)ds 


< 


< 


”■ rxi 

\ai{{Pi)t,s)\ ds 


^ rxi 


f3i{{qi)t, s) ds + 


^ rxi 


Wi{t, s) ds. 


Since the last summand in the above formula represents zero in 'M, this 
proves our conclusion for a, 6 G M. 

Now, we can proceed as in the proof of Thm. 29: 


f 


f 


= (6 - a) • 


< {b — a) 


f 
■f 


f [a + {b — a)a\ da 


g[a + {b — a)a] da 


□ 

We have just showed that Thm. 30 follows from Lem. 16 and the analo¬ 
gous property of the Riemann integral. In the same way, we can prove the 
following: 

Theorem 31 (Cauchy-Schwarz inequality). Let J be a non-infinitesimal 
interval of *M, and let f, g : J —*M he smooth functions. Let a, b ^ J, 
with a<b. Then fgf < (J^'' f) ■ (fl^g^). 
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9. Multiple integrals and Pubini’s theorem 


In previous sections, we have developed integral calculus with an interval 
as integration domain. In this section we will define the integral (and multi¬ 
ple integral) of a smooth function with a more general integration domain. 
The approach is similar to that used for the Peano-Jordan measure. 

We denote by J' the set of all the intervals of ’M (see Sec. 2). Thanks 
to the total order relation < on *M, every interval uniquely determines its 
endpoints. (We recall that this property is not true e.g. in SDG, where < 
is a pre-order.) Since the domain of every smooth function defined on an 
interval of *M can always be extended to a Fermat open subset (Lem. 8), we 
can therefore dehne the integral over an interval as follows: 


Definition 32. Let J be a non-infinitesimal interval, let / : J —>■ *M be a 
smooth function, and let I G be an interval contained in J. We define 




sup(/) 

inf{7) 


/• 


The following lemma permits to order the endpoints of pairwise disjoint 
intervals. 


Lemma 33. If I, J ^ ^ are disjoint intervals, then either sup{I) < inf(J) 
or sup(J) < inf (I). 

Proof. Set a := inf(/), b := sup(/), a := inf(J) and b := sup(J). By 
contradiction, we assume d < b and a <b. We distinguish four cases: a < d 
and b < b, a < d and b <b, d < a and b < b, d < a and b <b. In each case, 
we can find a point in / fl J. □ 


This lemma permits to integrate smooth functions over a finite pairwise 
disjoint union of intervals: 


Theorem 34. Let {Ii, ..., In} and { Ji, ..., Jm} be two finite pairwise dis¬ 
joint families of J' such that 

n m 

U = U 

i=l j=l 


Let J be 
function. 


a non-infinitesimal interval and let f : J 
Then 


n „ m ,. 

i=i 


/=://■ 


'U 


•M be a smooth 


(38) 
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Proof. Set Oi := inf(/j), 6j := sup(/j). Thanks to Lem. 33, we can always 
rearrange the indices i, j so that 

fll ^ ^ ^ Ofi ^ b^ri. 


Since the intervals are pairwise disjoint, an equality of the type bi = Oj+i 
may happen only if bi ^ li or Oj+i ^ li+i- If li U Ij+i is connected, we 
could replace li and Ij+i in the set {Ii ,..., /„} by li U li+i to get another 
hnite pairwise disjoint family of J^. In this case, because 


hi f + fli+i f - fah' f ~ hiUli+i f ’ 


we know that 


n „ n—1 „ 

hJ-J-U 


In other words, we can combine two “joinable” intervals to get a new set 
of hnite pairwise disjoint family, and keep the integral unchanged. After 
hnitely many such procedures, we reach a set of hnite pairwise disjoint family 
consisting of all the connected components of U, and by the invariance of 
integrals, we get 


U'-Pl 


/ I / is a connected component of U 


We arrive at the same result by considering YlJLi fj f- 

In the following dehnition, we identify the general domains for integrals 
of higher dimensions. 

Definition 35. A box in is a set of the form B = \hi=ih-, where 
are intervals. We denote by the set of all the boxes in •R'^. An 
elementary set U € S' is a, set of the form ljr=i where Bi G are boxes. 


Theorem. S is a Boolean algebra with respect to union, intersection and 
set difference. 

Proof. It is enough to show that S is closed under these operations. The 
closure of S with respect to union is clear. HU = UILi ^ ~ Uk=iCk, 

where Bi = and Ck = 0^=1 Jkj are boxes, then 


r d 

Ur\V = \JlW^Iijr\Jkj\i = l,...,n,j = l,...,m 
[i=i 

Since fj n Jkj G (including the empty interval), U UV G S. To prove 
the closure of S with respect to set difference, we hrst introduce the notion 
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of complement := \ A for any subset A C •R'^. Using the same 

notations, since 




( d (l-l d 'I 

n = U m'M X X J] 'M I / = 1,..., s = ± U 

j=l ) U=1 i=l+l ) 


where = {x G ’M | x > y Vy € Jki} and {J^i) = {x G ’M | x < 

y 'iy ^ Jki}- we know that G (f, and hence 

U\V = U es. □ 


In our framework, the following result corresponds to Fubini’s theorem. 

Theorem 36. Let B = OiLi ^ ® ^ ^ 

a Fermat open set such that B C U, and let f : U —>■ *M be a smooth 
function. Let (ui,..., be a permutation of {1,..., d}. Then 



We denote this integral by f- 

Proof. We prove only for the case d = 2, since the general case is similar. Let 
a, 6 , a, b be the endpoints of Ii and I 2 , respectively. Set F{t) := f{t, y) dy 
and G{t) := ^ /(x, y) dx^ dy for each t G Ii. Both F and G are smooth 

functions such that G(a) = 0 and 

G'(t) = y^ ^ (/ = = F(t) vtG/i. 

Therefore, G = F and G{b) = (^fl^f{x,y)dx^ dy = /V(t) df = 
fa {id fi^^y) dy) dx, which is our conclusion. □ 

Similarly to Thm. 34, we can finally prove the following: 

Theorem 37. Let {Bi,..., Bn} and {Ci,... ,Cm} be two finite pairwise 
disjoint families of such that 

n m 

UBi=UCj=:t/. 

i=l j = l 

Let V be a Fermat open subset of containing U, and let f : V —*M 
be a smooth function. Then 
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10. Conclusions 


The main aim of the present work is to lay the foundation of the integral 
calculus for smooth functions defined on the ring of Fermat reals. In the 
Introduction 1, we have already summarized the major achievements of this 
theory compared to other non-Archimedean theories. We can outline them 
by saying: better order properties and easier to deal with infinite dimensional 
integral operators. 

The present work confirms that the ring of Fermat reals does not repre¬ 
sent a new foundation of the calculus. On the other hand, the dehnition of 
integral as primitive (Def. 14) and the corresponding simple “global proofs” 
(see e.g. Thm. 17) can suggest ideas for a new approach to a subset of the 
classical integral calculus. On these bases, we can also think of the possibil¬ 
ity to approach the integral calculus of smooth functions from an axiomatic 
point of view. 

Finally, our work shows that all the classical instruments to deal with in¬ 
tegrals of smooth functions with infinitesimal parameters are now available. 
This is important for the applications in Physics. 

A suitable notion of convergence for sequences of smooth functions, the 
exchange between limit and integral, and a more general class of integration 
domains remain open. We plan to explore them in a future work. 
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